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Abstract 

It is shown that MSSM with first two generations of squarks and sleptons much heavier 
than the third one naturally predicts the maximal stop mixing as a consequence of the RG 
evolution, with vanishing (or small) trilinear coupling at the high scale. The Higgs boson is 
generically heavy, in the vicinity of 125 GeV. In this inverted hierarchy scenario, motivated 
by the supersymmetric FCNC problem and models for fermion masses based on horizontal 
symmetries, the mass of the lightest stop is 0(0.5) TeV and of gluino - 0(2 — 3) TeV. The LSP 
can be either higgsino or bino or a mixture of both and it can be a good dark matter candidate. 
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1 Introduction 



The issue of naturalness of the Higgs potential has for long been a driving principle for the- 
oretical ideas going beyond the Standard Model. This is a qualitative and only a theoretical 
argument. However, if one abandons it, there is no reason to expect new physics at the LHC. If 
one takes it too literally, new particles should have already been observed at the LHC or even 
at LEP. 

In supersymmetric (SUSY) models, MSSM in particular, it has been appreciated from the 
beginning of the supersymmetric phenomenology that only part of the superpartner spectrum is 
concerned by the naturalness argument, namely higgsinos, the third generation sfermion masses 
and gluinos (and to a much lesser extent other gauginos). Those are the particles whose masses 
enter the Higgs potential either at the tree level (higgsinos) or via loop effects, enhanced by 
large couplings. Naturalness argument suggests that those particles are rather light but says 
nothing about the masses of the first and the second generation sfermions. 

On the other hand, it has very early been observed that the supersymmetric FCNC and 
CP violation problems can be substantially eased the heavier the first-two-generation sfermions 
are. Thus the concept of naturalness and the FCNC supersymmetric problem lead together to 
the expectation of a split sfermion spectrum. Furthermore, in fermion mass models based on 
horizontal symmetries, the sfermion mass spectrum has been linked to the hierarchical fermion 
masses. The predicted pattern is the so-called inverted hierarchy (IH) of the sfermion masses. 

In this paper we investigate the predictions for the lightest Higgs boson mass in the MSSM 
with inverted hierarchy of sfermion masses. Due to 2-loop effects, heavy first two generations 
play significant role in reaching the Higgs mass in the region of 125 GeV which seems to be 
favoured by the recent ATLAS [lj and CMS [2] data. This is because the maximal stop mixing 
is obtained from the RG evolution effects, with initial A < mo (3) and even with A = 0, with 
mo (3) being the scale of the 3rd generation sfermion masses. The lightest stop is predicted to 
be in the 500-1000 GeV mass range and gluinos are 2-3 TeV heavy. The model is less fine-tuned 
than CMSSM. 

Various phenomenological aspects of the inverted hierarchy of the sfermion masses and its 
impact on the Higgs mass have also been discussed in a recent paper |3j. Our findings are in a 
qualitative agreement with the results of ref. Our predictions for the Higgs boson mass are 
often a couple of GeV higher and, as we discuss in more detail at the end of Section [3j this can 
be traced to the 1-loop [3] versus 2-loop (in our case) calculation of the Higgs mass and to the 
not fully overlapping range of the investigated parameter space. 

The paper is organized as follows. In section (J2J) we discuss the implications of the low 
energy MSSM spectrum on the Higgs boson mass. A crucial role of the stop mixing in reaching 
large values of the Higgs boson mass is emphasized. In section ([3j) we perform a detailed study 
of the predictions for the Higgs boson mass in the IH scenario. We also review there some of 
the theoretical ideas leading to the IH of sfermion masses. In section [4] other phenomenological 
constraints and implications for the LHC are discussed. Several benchmark points are presented. 
Our conclusions are presented in section |5j 
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2 The Higgs boson mass in the electroweak scale MSSM 



It is well known that at the tree level the mass of the lightest Higgs boson in the MSSM is 
bounded from above by Mz\ cos 2/3 1. However, the tree-level upper bound on the Higgs boson 
mass is uplifted at the quantum level when the effects of (soft) SUSY breaking are taken into 
account. The magnitude of the loop corrections depends mainly on the properties of the stop 
sector. In the decoupling limit, 3> M z where A denotes the CP-odd scalar, the Higgs boson 
mass corrected by the dominant one-loop contribution is given by [4]: 
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renormalization scheme) and X t = A t — /x/tan/3 with A t being SUSY breaking top trilinear 
coupling at Msusy- It follows from the above formula that sizable corrections to the Higgs 
boson mass can be obtained when the stops are substantially heavier than the top quark. Since 
the contribution to the Higgs boson mass from stop mixing may be significant, the precise 
values of the stop masses required to obtain a given value of the Higgs boson mass are quite 
sensitive to the ratio X 4 /Msusy- It follows from eq. ([I]) that the maximal contribution from 
stop mixing is obtained for (XjI/Msusy = \/6. 

Equation ([I]) was derived under the assumption of a relatively small mass splitting between 
the stops. The generalization of this formula to the case of a large splitting can be found e.g. in 
Ref. [Ijj. One can then show that, for fixed Msusy, the value of \X t \ /Msusy giving the maximal 
correction from stop mixing to the Higgs mass, as well as the value of this maximal correction 
increase with the splitting between the two eigenstates of the running stop mass matrix. 

There are theoretical uncertainties in the prediction of the Higgs boson mass. They originate 
from unknown higher order corrections and from the limited experimental knowledge of the SM 
parameters, mainly the top mass and the strong gauge coupling constant. The estimated 
uncertainty is about 3 GeV [Bj Q Taking this into account the Higgs boson masses in the range 
122-128 GeV are consistent with the 125 GeV Higgs interpretation of the excesses observed by 
the LHC experiments. In the present paper we use mf ole = 173.3 GeV and a s (M z ) = 0.1187. 

In our numerical studies we use a modified version of SOFTSUSY v3.2.4 [8] which employs 
the two-loop formulae for the Higgs boson mass [9]. The original version of SOFTSUSY refuses 
to calculate the spectrum if there are any negative running {DR) squared masses at the Mz 
scale, even if there are no tachyons at the Msusy scale where the Higgs potential is minimized. 
The reason for this is that in SOFTSUSY the running parameters are used to compute the 
SUSY corrections to the SM parameters (such as gauge or Yukawa couplings) at the M z scale. 
Such a procedure leads to excluding some of the parameter space which is perfectly viable 
from the theoretical point of view. Since this part of the parameter space is important for the 
studies in the present paper we have modified SOFTSUSY in such a way that the pole masses 
are used in the calculation of the SUSY corrections to the SM parameters and the tachyons are 
signalized only if there are any negative squared running masses at the Msusy scale. 



1 The dominant 3-loop contribution to the Higgs mass has been calculated in [6 . In the case of CMSSM it 
is positive and of the order of 1-3 GeV with precise value depending on the soft terms [7]. 
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Figure 1: The Higgs boson mass versus Msusy for various values of X t /M$ and tan/3 = 10 
(top) or tan/3 = 50 (bottom). The solid (dashed) lines correspond to the negative (positive) 
values of X t . All the other MSSM parameters (defined at the Msusy scale) are fixed to be 2 
TeV except for M 2 = fi = rriA = 1 TeV, Mi = (to ensure that a neutralino is the LSP), 
Af, = A T = 0. Analogous plots for negative \i — — 1 TeV are very similar. 



In Figure [T] we plot the Higgs boson mass versus Msusy for various values of X t /Ms, where 
Ms = yJuiQirnu and are the running soft masses of the third generation left-handed 

squark doublet (right-handed up-type squark) at the M S usy scalej^] The data used in this 
plot were obtained by scanning Ms = vhq = mjj up to 5 TeV while keeping fixed all other 

2 The parameter Ms is defined using the soft masses while Msusy introduced in Eq. ([lj is defined using the 
eigenvalues of the running stop mass matrix. The values of Ms and Msusy ar e quite similar in most part of 
parameter space. The relative difference between Ms and Msusy increases as Mg decreases and in some cases 
may reach 10%. In some figures we use Ms for technical reasons. In SOFTSUSY it is possible to fix the value 
of Ms while Msusy is obtained as a result of an iteration procedure. 
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Figure 2: The Higgs boson mass versus X t /M$ for various values of the ratio mq/mu keeping 
fixed Ms = 1 TeV assuming tan/3 = 10. All other MSSM parameters at Msusy are the same 
as in the case of Figure [[} 



MSSM parameters defined at the Msusy scale. It is clear from this figure that the scenario 
with vanishing or very small stop mixing is incompatible with the Higgs boson mass of about 
125 GeV (suggested by recent ATLAS [I] and CMS [2] results). The Higgs boson mass of 120 
GeV may be obtained in the absence of the mixing in the stop sector but only for very heavy 
stops. On the other hand, when the contribution to the Higgs boson mass from the stop mixing 
is large, stop masses of about 1 TeV may be consistent with the Higgs mass of 125 GeV. The 
largest Higgs mass is obtained for \X t /Ms\ between 2 and 2.5, in agreement with the leading 
one-loop formula (JIJ. Notice also that even though the expression (JTJ does not depend on the 
sign of X t , for the positive X t the Higgs mass is enhanced more than for the negative X t . This 
asymmetry arises due to SUSY threshold effects on the top Yukawa coupling which depend on 
the product of the gluino mass and A t [TO] . 

The dependence of the Higgs mass on the stop masses splitting is shown in Figure [2] The 
Higgs mass versus X t /M s is plotted there for various values of the ratio mq/mjj keeping fixed 
the value of Ms = 1 TeV. While for mq = mjj the local maxima of the Higgs mass oc- 
cur at | X t /Ms | m ax ~ 2.2 the positions of these maxima move towards the larger values of 
\X t /Ms\ ma , x exceeding 3 for mq/mu ^ 5. Moreover, the value of the Higgs mass corresponding 
to (Vt/Mslmax in the case of large stop mass splitting may be larger by several GeV as compared 
to the unsplit case. Figure [2] shows that the bigger stop mass splitting is the bigger maximal 
Higgs mass may be obtained. It is also interesting to note that for a given \X t /Ms\ < 2 the 
mass splitting between the stops tends to suppress the Higgs mass. 

For large values of tan /3 for which the bottom and tau Yukawa couplings are of the same 
order as the top Yukawa coupling the loop corrections from the sbottom and stau sectors may 
also be significant. These corrections tend to reduce the Higgs boson mass and become quite 
large for \fi\ ^> and/or ^> m f , especially when /zM 3 < since this typically leads to the 
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enhancement of the bottom Yukawa coupling due to SUSY threshold corrections important for 
large tan /3. These effects have been recently discussed in [TT]. However, such reduction of the 
Higgs mass at large tan (3 is not present in Figure [I] (also for Figure [2] the (not shown) plots for 
large tan /3 are very similar to those for tan /3 = 10) which is a consequence of the fact that |//| 
is chosen there to be smaller than the sbottom and the stau masses. For the same reason the 
results for negative /x are very similar to those for positive \i presented in these Figures. As a 
matter of fact, for our choice of the parameters the Higgs mass for tan (3 = 50 is slightly larger 
than for tan = 10. This can be attributed to the tree-level contribution which grows with 
tan (3. 

3 The Higgs boson mass in IH Model 

As discussed in the previous section, Msusy > 1 TeV and a big contribution from the stop 
mixing to the Higgs mass is necessary to obtain rrth > 122 GeV. The soft trilinear term A t at the 
EW scale should not be much different from its optimal value corresponding to |X t /MgusY|max 
which maximizes the Higgs mass (for fixed other parameters). Typically, such optimal value 
of | A t | is 2 to 3 times bigger than Msusy- In UV models, the values of Msusy and the stop 
mixing at the EW scale depend on the pattern of soft terms fixed at a (usually high) scale of 
supersymmetry breaking mediation and on the RG evolution. In this paper we are interested 
in the prediction for the lightest Higgs boson mass in the scenario with the inverted hierarchy 
of squark masses, motivated by the naturalness arguments and by the flavour models based 
on horizontal symmetries. We shall not refer to any specific model. The Higgs boson mass is 
calculated for a given set of the boundary conditions for the soft terms at the GUT scale. 

Since the large stop mixing at the EW scale is crucial for the successful prediction for the 
Higgs boson mass, before going to the IH scenario, it is useful to discuss when large values of 
| A t /M SU sY I at the EW scale can be obtained. Due to the RGE running, the low scale value of 
At is a linear combination of the high energy scale gaugino masses and trilinear terms. In the 
case of the universal high energy scale boundary conditions for these parameters (M1/2 and Aq, 
respectively), the EW scale value of A t is approximately given by 

A t w -1.6Mi /2 + 0.354, • (2) 

The precise values of the coefficients in the above equation depend on the gauge and Yukawa 
couplings, so implicitly (among other parameters) on the values of the top mass and tan /3. The 
coefficient in front of Aq is relatively small because the top quark is heavy. It goes to zero when 
the top quark mass approaches its infra-red quasi-fixed point value [121 f!3] . 

Equation (|2]) is valid at the 2-loop RGE level. For the sake of definiteness we give in this 
paper all the coefficients obtained from the RG running (unless stated otherwise) computed 
at the scale Q = 1.5 TeV and the gauge and Yukawa coupling at the GUT scale consistent 
with radiative electroweak symmetry breaking (REWSB) for tan/3 = 10, Mi/ 2 = 700 GeV, 
A = —3 TeV, the scalar mass of the third generation m (3) = 3 TeV and that of the first and 
second generation m (l, 2) = 10 TeV (A t in eq. Q does not depend on the latter). In most of 
the paper, unless otherwise stated, we also take universal Higgs soft masses at the GUT scale 
m Q (H u ) = m (H d ) = m (3). 
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Figure 3: Contours of the Higgs boson mass (black dashed line), the lighter stop mass (solid 
green line) for /i > 0, tan/3 = 10, iWi/ 2 = 1 TeV and A = —2 TeV. The yellow "tachyonic 
stop" and the grey "no REWSB" (/i 2 < 0) regions are excluded. In the dark green region the 
relic density of neutralinos gives Odm^ 2 < 0.1288 [16J. 

The values of niq and rrijj at the EW scale in terms of the GUT scale values of the soft 
terms are given by: 

m 2 Q « 3.1M? /2 + 0.1A M 1/2 - 0.04A 2 + 0.65m (3) 2 - 0.03m (l, 2) 2 , (3) 
ml « 2.3M 2 /2 + 0.2A M 1/2 - 0.07A 2 + 0.35m (3) 2 - 0.02m (l, 2) 2 . (4) 

The dependence of the 3rd generation squark masses at the EW scale on the m (l, 2) is a 2-loop 
effect and it is negligible for small values of m (l,2). However, as was pointed out in [TJ], for 
very large scalar masses of the first two generations it may lead to tachyonic stops due to the 
(small) negative coefficient in front of mo(l,2) 2 in the expressions Q-Q. 

An inspection of Eqs. ([2])-(j4]) shows that, as long as the 2-loop effect is negligible, e.g. in 
CMSSM, "maximal" stop mixing can be obtained only for sufficiently large values of A . The 
needed value of A depends on the relative magnitude of Mi/ 2 and mo- For m ^> Mi/ 2 , 
Ao/mo ~ ±2 is required, with negative sign more effective for not too large a ratio mo/M]_/2- 
For m <C My2 one needs A /M 1 / 2 ~ —3.5 and for Mu 2 ~ m that ratio has to be about —4. 
These conclusions agree with a recent study [15J. The reason a rather large \A Q \ is needed is 
that the RG evolution of the dividend and of the divisor in the ratio Af/(mQmu) is correlated 
(leading to A 2 /M| USY < 1 for small \Aq\). In consequence, a 125 GeV Higgs boson requires a 
heavy spectrum and significant cancellations in the Higgs potential. 

In the IH scenario, with m (l,2) 3> ^o(3), the RG evolution of A t can be decoupled 
from the evolution of stop masses because the former does not depend on mo(l,2). One can 
enhance A t by gluino contribution to the RG running, without enhancing the stop masses, due 
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Figure 4: Slices of the part of parameter space shown in Figure [3] for several fixed values of 

m (3) given in TeV on the plots. The solid (dashed) lines correspond to the Higgs mass (the 

lighter stop mass, M SU sy ot X t /M SU sY, respectively). 
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Figure 5: The same as in Figure [3] but for M 1 / 2 = 1.5 TeV and A = 0. 



to the negative 2-loop contribution to the stop masses from the 1st and 2nd generations. In 
consequence no large initial values of Aq are needed for obtaining the Higgs mass in the 125 
GeV range, with the lighter stop mass in the range 500-1000 GeV and the gluino mass 2-3 TeV. 
Thus, Msusy ~ 1 TeV and "maximal" mixing follow quite naturally from the RG evolution. 
Moreover, the cancellations in the Higgs potential are significantly smaller than e.g. in the 
CMSSM focus point. 

An example of the predictions for the Higgs boson mass in the IH scenario, for tan = 10, 
Mi/2 = 1 TeV and A = —2 TeV is shown in Figure [3j Very similar predictions are obtained 
for a range of M\m between 0.5 and 2 TeV and A between zero and about —4 TeV, if larger 
Mi/2 is taken with smaller \A \. 

It can be seen from Figure [3] that IH scenario predicts large values of the Higgs boson 
mass. In order to better illustrate the dependence of the Higgs mass across parameter space, 
particularly near the tachyonic region, we present in Figure [4] several slices of the plot from 
Figure [3] for three fixed values of m (3). For a given value of m (3), the stop masses, as well 
as Msusy; decrease with increasing mo(l,2). In consequence, the logarithmic correction to the 
Higgs mass decreases while the one from stop mixing increases. Since the latter correction 
increases polynomially with \X t \/Ms\jsY the Higgs mass initially increases reaching a maximum 
for a value of mo(l,2) corresponding to |Xj|/Msusy ~ 2. Increasing m (l,2) further decreases 
rrih because the stop-mixing correction is close to its maximal value and its increase cannot 
compensate the decrease of the logarithmic correction. Moreover, for large enough |X t |/MgusY 
the stop-mixing correction also starts to decrease resulting in a rapid decrease of the Higgs 
massJ3 

3 A sharp cut-off of some curves is partly due to the fact that SOFTSUSY for large mo (3) is not able to reach 
very small values of Msusy before the tachyonic stop is signalized. 
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Figure 6: Slices of the part of parameter space shown in Figure [5] for several fixed values of 
mo (3). The solid (dashed) lines correspond to the Higgs mass (the lighter stop mass). 



Non-trivial constraints on the IH scenario follow from the requirement of proper REWSB. 
The origin of these constraints can be understood by inspecting a dependence of /i 2 on the soft 
terms defined at the GUT scale which is approximately given by: 



fi 2 » -m 2 Hu « 1.3M 2 /2 + O.LA 2 - 0.35M 1/2 A) - 0.01m (3) 2 - 0.006m (l, 2f 



(5) 



Notice first the smallness of the coefficient in front of mo(3) 2 . Even though it is negative 
in the above formula we should emphasize that its sign depends on the scale at which the 
coefficients are extracted and to some extent on the boundary conditions for the soft terms 
(which influence the GUT scale values of gauge and Yukawa couplings). It is important to 
note that the coefficient in front of mo(3) 2 becomes more and more negative as the Higgs 
potential minimization scale, Msusy, increases. On the other hand, the very small coefficient 
in front of m (l,2) 2 is always negative due to the specific structure of the two-loop RGEs. 
Even though this coefficient is very small, large values of mo(l,2) may drive /i 2 negative. In 
order to protect proper REWSB the positive contribution from Mi/ 2 and A has to be large 
enough to ensure that by increasing mo(l, 2) the stops become light enough to realize "maximal 
mixing" scenario before /z 2 becomes negative. This effect is seen in Figure [3j The "maximal 
mixing" can be realized only for mo (3) below some critical value which in this case is smaller 
than 5 TeV. For larger values of m (3), fi 2 becomes negative before m (l,2) reaches the value 
corresponding to (|X t | /Msusy) max which gives the maximal stop-mixing correction to the Higgs 
mass. Constraints from REWSB are also the reason for the lack of maximum in a dependence 
of the Higgs mass on mo(l, 2) for a large values of mo(3), as seen in Figure |4j 

It is also evident from Figure [3] that the largest Higgs masses correspond to the largest 
values of m (l, 2). Therefore, a heavier Higgs prefers the regions of parameter space where the 
first and second generation sfermions are decoupled which makes plausible the solution to the 
FCNC problem. This also can be seen from the opposite perspective. If the solution of the 
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Figure 7: The same as in Figure [3] but for tan/3 = 50 and m# d = 1.6mo(3). The region below 
the purple line is excluded by BR(B S — > // + at 95% C.L. The orange region is excluded 
because it predicts a tachyonic stau. 



flavor problem relies on the decoupling of the first two generations of sfermions the Higgs mass 
is expected to be large. In particular, for the parameters used in Figure [3] the requirement 
of mo(l,2) > 15 TeV implies that the lightest Higgs mass has to be not smaller than about 
123 GeV. More generally, we found that if the lighter stop mass is at least 0(500) GeV and 
mo(l,2) > 15 TeV the lightest Higgs is necessarily heavier than about 122 GeV. This lower 
bound on the Higgs mass becomes even more stringent for larger values of mo(l, 2). 

Figure [5] demonstrates that in the IH scenario the stop- mixing correction to the Higgs boson 
mass can be maximized even if Aq = 0. In such a case A t at the EW scale is generated entirely 
radiatively by the gaugino masses. This requires gaugino masses somewhat heavier than in the 
case of non-vanishing A-terms. We found that for A = 0, = 122 (125) GeV can be obtained 
for Mi/2 > 900 (1500) GeV. Notice that in this example the Higgs mass of 124 GeV can be 
obtained even for very small values of mo(3)|^] This is because the value of M\/i = 1.5 TeV 
is large enough to generate stop masses radiatively. Smaller values of mo (3) require smaller 
?7io(l, 2) to avoid tachyonic stops. Also Figure [5] shows that the requirement of maximal mixing 
and of proper REWSB may put an upper bound on m (3). Dependence of the Higgs and the 
lighter stop mass on mo(l, 2) for Aq = is qualitatively similar to the one with non- vanishing 
A-terms, as seen in Figure [6j 

The IH scenario can be also realized for large values of tan (3. In Figure [7] the plot analogous 
to that from Figure [3]but for tan /3 = 50 is presented. In addition, the splitting between the soft 
Higgs masses at the GUT scale, m (H u ) = m (3), m (H d ) = 1.6m (3), is introduced in order to 

4 One cannot, however, take mo(3) to be arbitrarily small because otherwise the stau would be lighter than 
the lightest neutralino. 
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Figure 8: Slices of the part of parameter space shown in Figure [7] for several fixed values of 
m (3). The solid (dashed) lines correspond to the Higgs mass (the lighter stop mass). 



reduce the impact of the BR(£> S — > /i~V~) constraint on the parameter space. Even though a 
dependence of the Higgs and the lighter stop masses on m (l, 2) for large tan/3 is qualitatively 
similar to that for tan/3 = 10, as seen in Figure |8j the upper bounds on m (3) and mo(l,2) 
set by the condition of positive fi 2 are more stringent for tan /3 = 50. There are two reasons 
for that. First: the positive contribution to fi 2 from Mu^ and Aq is smaller for larger values of 
tan/3. Second: mo(ffd) (which is now larger than m (3)) gives larger negative contribution to 
fi 2 . Since the tau Yukawa coupling is now of the same order as the top Yukawa it gives large 
negative contribution to the stau masses and some part of parameter space at smaller values of 
m (3) is excluded because it leads to a tachyonic stauj^] However, it can be seen that in spite of 
these constraints the Higgs boson mass may reach even slightly larger values for tan (3 = 50 as 
compared to the tan /3 = 10 case. Moreover, gaugino masses required to obtain a given value of 
the Higgs mass are somewhat smaller for tan /3 = 50. In particular, for A = the Higgs boson 
mass of 122 (125) GeV can be reached if M 1/2 > 800 (1300) GeV. 

The examples presented above were designed in such a way that the Higgs mass of 125 GeV 
is obtained with the smallest possible values of M 1 / 2 and \A \. We should, however, stress that 
in the IH model larger Higgs masses can be reached if Mi/ 2 and/or \A \ are chosen to be larger. 
For example, for Mi/ 2 = 2.5 TeV and A = the Higgs mass up to about 127 GeV can be 
obtained. 

One of the reasons why the "maximal mixing" generated by IH is interesting is that it 
requires less fine-tuning than e.g. CMSSM to obtain the same Higgs boson mass. In order to 



5 As a matter of fact, for tan/3 = 10 the tachyonic stau is also present in some part of parameter space but 
for much smaller values of mo (3) (outside of the range of the plot in Figure [3]). 
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quantify this we define the fine-tuning parameter A = max{A a }, where j^] 

d In rrih 



A, 



d In a 



(6) 



The index a stands for any soft term or /x. For a given point in parameter space proper REWSB 
requires cancellation between the parameters with the precision of order 1/A. As explained 
before, in IH scenario with the "maximal stop mixing" the bigger the value of M1/2 is, the 
smaller value of Aq is needed to obtain a given value of the Higgs mass. In other words, the 
same Higgs boson mass can be obtained for various values of the ratio Ao/M^-H The level of 
fine-tuning strongly depends on the relative values of Aq and M\i%. It can be inferred from the 
coefficients in eq. (|5]) that from the naturalness point of view neither large nor vanishing Aq is 
the best choice. We found that for tan /3 = 10 and A = the Higgs mass of 125 GeV requires 
A ~ 0(250) but for e.g. A Q = -2 TeV and M 1/2 = 1 TeV the Higgs mass of 125 GeV can be 
reached with A ~ 0(150). Similarly, for tan/3 = 50 the 125 GeV is reached with A ~ O(150). 

It is also interesting to note the existence of a generalized "focus point" region. We recall 
that in the CMSSM the "focus point" refers to the region of parameter space where ttiq 3> M1/2 
[T9| 120] . The focus point in the CMSSM gained a lot of interest since heavy stops required to 
lift the Higgs mass above the LEP bound [21] can be reconciled with small contribution of m 
to the Higgs potential and keeping /i 2 close to the EW scale. Such small values of /i 2 result 
from the cancellation of the gaugino contribution to m 2 Hu with that of scalars (see Eq. (JsJ) ) . In 
the focus point region of the CMSSM the LEP bound can be satisfied without introducing large 
fine-tuning of the EW scale. However, in the context of the 125 GeV Higgs the focus point 
region of the CMSSM is much less attractive because in this part of parameter space the stop 
mixing is strongly suppressed, unless Aq is large [22} [23]. In the IH parameter range, by taking 
mo(l,2) >> mo(3) > Mi/ 2 , the virtues of the focus point remain but in addition one obtains 
large stop mixing and the region of 125 GeV for the Higgs mass is easily reachable. 

Our results are in a qualitative agreement with the results of ref. [3]. Somewhat smaller 
values of m h reported in ref. [3] are due to the fact that in ref. [3] the computation of the 
Higgs boson mass is performed at one-loop level using ISASUGRA program from the ISAJET 
package [21], while we use the routines adopted in SOFTSUSY [8] which take into account all 
relevant two-loop corrections, n For example for the benchmark points listed in Table 1 of ref. 
[3] we find values of m h larger by 0.5-2 GeV. Secondly, in ref. [3] the region of light (below 1 
TeV) average mass fhq(3) = (m^ + mj 2 + mgj /3 is emphasized as the one that satisfies simple 



6 In the case of large loop corrections to the quartic Higgs coupling, it is more appropriate to refer to 
rather than Mz when defining the measure of fine-tuning, see e.g. Refs. [T7l I18j . 

7 Notice that a given value of the Higgs mass can be obtained only for the values of the ratio of Aq jM-y / 2 in 
some finite range whose size is controlled by the value of Mm- 

8 There is also another difference between ISASUGRA and SOFTSUSY. In ISASUGRA the heavy sfermions 
of the first two generations are decoupled from the RGEs at a scale equal to the heavy-sfermion mass, while 
SOFTSUSY decouples all the SUSY particles at a common scale Q = M z . In SOFTSUSY only leading 
logarithms of the ratios between the heavy-sfermion masses and Mz are included in the threshold corrections 
to the gauge and Yukawa couplings. On the other hand, the decoupling procedure adopted in ISASUGRA 
effectively accounts for the resummation of the logarithms. Simple estimate of the contribution from the 
non-leading logarithms, which is missing in SOFTSUSY, gives few percent for the sfermion masses of order 
O(lOTeV). 
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naturalness constraints. Our values of fhq(3) are typically above 1 TeV (because of the masses 
of the heavier stop and the lighter sbottom) and for the average rhq(3) in the range 1-1.5 TeV 
the upper bound on rrih reported in ref. [3] is 124 GeV, quite close to our result of 125-126 
GeV. We also notice that the fine tuning in our parameter range, although with heavier third 
generation spectrum, is similar to that for the benchmark points of Table 1 of ref. [3] , because 
of our generalized focus point. 

3.1 Inverted hierarchy models: Abelian flavor symmetries 

So far, we have discussed the predictions for the Higgs boson mass in the IH scenario, with 
m (l, 2) >> m (3), without any reference to potential models of such hierarchy. We would like 
now to recall briefly its link to models of fermion masses based on horizontal symmetries. 

Inverted hierarchy was proposed some time ago [25j [26] as a way to ease the FCNC and CP 
constraints in supersymmetric models. Early ideas did invoke horizontal non-Abelian symme- 
tries for explaining fermion mass hierarchies, like U(2) [25], under which first two generations 
transform as a doublet, whereas the third generation is a singlet. Whereas U(2) models do 
explain the difference between the first two generations and the third one and therefore can 
accommodate a hierarchy between the first two and the third generation of scalars, they do 
not actually predict it. To our knowledge, the first class of models in which the inverted hier- 
archy was really predicted [27] are supersymmetric generalizations of Abelian flavor models of 
the Froggatt-Nielsen type [28J. These models contain an additional Abelian gauge symmetry 
U(l)x under which the three fermion generations have different charges (therefore the name 
horizontal or flavor symmetry), spontaneously broken at a high energy scale by the vev of (at 
least) one scalar field $, such that e = (<&)/M << 1, where M is the Planck scale or more 
generically the scale where Yukawa couplings are generated. Quark mass matrices for example, 
in such models are given, order of magnitude wise, by 

ij ' ij ' 

where q^ (ui,di,h u ,hd) denote the U(l)x charges of the left-handed quarks (right-handed up- 
quarks, right-handed down-quarks, H u and Hd, respectively). 

A successful fit of the experimental data requires larger charges for the lighter generations 

qi > <?2 > <?3 , ui > u 2 > u 3 , di > d 2 > d 3 , (8) 

one simple example being (see e.g. |29|) 

(qi, Mi, ei) = 3 , (q 2 , u 2 , e 2 ) = 2 , (q 3 , u 3 , e 3 ) = , (d h h) = 1 , (d 2 , l 2 ) = , (d 3 , l 3 ) = . (9) 

Scalar soft masses in Abelian flavor models are typically of the form (only the diagonal in 
flavour entries are relevant for our present discussion) 

m\ = Qi(D) + (mf) 2 , (10) 

where Qi(D) are D-term contributions for the scalar of charge Qi, whereas {mf) 2 are F-term 
contributions. D-term contributions were argued to be naturally generated (at least) in effective 
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string models [30], to be positive and, in certain circumstances, to be dominant over the F-term 
contributions. It is then clear from (j8|, ^ and (10) that the hierarchy of the masses of scalars 



is predicted to be inverted with respect to the hierarchy of fermion masses. 

Whereas Abelian models naturally predict the inverted hierarchy, they do not generically 
predict approximate degeneracy among the first two generations, unlike their non-Abelian 
cousins. This leads to possible tension with FCNC constraints, which have to be analyzed 
in some details (for previous works see e.g. [27j EH])- We would like also to point out that re- 
cently there were various other explicit realizations of the inverted hierarchy, based on geometric 
localization and non-Abelian family models |31j . 

From the perspective of our present paper, the inverted hierarchy models do generically 
predict also a splitting between the first-two-generation sfermions as well as intragenerational 
splitting 

™ Qi ^ m Ui ^ m Di . (11) 

A relevant question is therefore what changes are expected in the results presented so far by 
relaxing the hypothesis of degenerate first two generations and of rriQ i = mjj. — m^. Since 
the first two generations are very heavy, we could expect large RGE effects. The RGEs of all 
scalar soft masses and in particular of the third generation of squarks and of the Higgs scalars 
depend at 1-loop level on the combination (see e.g. [32] ) 

3 

S = Tr(Ym 2 ) = m 2 Hu - m 2 Hd + ^[mj. - 1m\ + m 2 D . - m\. + m 2 E .) , (12) 

i=i 

where the trace is over the whole spectrum of MSSM states and which, under our assumption 
until now, is zero at high-energy. Interestingly enough, in Abelian flavor models with D-term 
dominance of the type put forward in this section, the quantity S is equal to 

S = Tr(FQ) (D) . (13) 

However, Tr(FQ) has to vanish (or to be very small) for phenomenological reasons, as argued 
in various papers [33] and therefore our conclusions concerning the running of soft terms and 
the fine-tuning remain unchanged. In particular, for the charge assignment ([9]) 5* = and the 
prediction for the Higgs boson mass are almost identical to those presented in the previous 
plots with degenerate first two generations. 

If the U(l)x charges of fermions are not chosen in such a way that S vanishes then S is 
generically of order (D) i.e. of the same order as the first two generation soft masses. Since S 
enters RGEs at the one-loop level its effect on the RG running of soft scalar masses (at least for 
non-colored scalars) would generically be much larger than the two-loop effect from heavy first 
two generations. The contribution from S to the EW scale soft scalar masses is determined by 
the hypercharge assignment and is approximately given by: 

m 2 f = -OMYfS, (14) 

where Yj is the hypercharge of fermion /. In particular, S contributes to the EW scale value 
of rn 2 Hu w —0.0255 (compare with Eq. (5)) so for the values of y/S ~ O(lOTeV) it is the 
dominant contribution implying that proper REWSB is possible only if S is positive. For 
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S > 0, S gives positive contribution to mfj ps 0.035S 1 and compensates the negative two-loop 
effect (see coefficients in Eq. Q) which is crucial for obtaining the maximal stop mixing. The 
contribution to mq ~ — 0.008S 1 is negative but relatively small. The largest contribution from 
S is received by m 2 E ps —0.055*. Since this contribution is negative one can expect that large 
values of S typically result in tachyonic staus. Therefore, we conclude that in IH models the 
maximal stop mixing consistent with REWSB is possible only if 15*1 is smaller than the D-term 
contribution to the first and second generation squarks. 

Let us also briefly comment on the effect of the mass splitting within the third generation 
of sfermions. Such splitting may be present due to arbitrary 0(1) coefficients of the diagonal 



F-term in eq. (10). The individual contributions from the scalars to /i 2 is approximately given 
by: 

/i 2 « -0.6m (# u ) 2 + 0.35mo(Q3) 2 + 0.25mo(£/ 3 ) 2 

-0.025 m (H d ) 2 + 0.025 m (D 3 ) 2 - 0.025 m (L 3 ) 2 + 0.025 m (E 3 ) 2 . (15) 

For universal soft scalar masses, the overall contribution from the scalars to fi 2 is quite insen- 
sitive to the scalar mass mainly because the contributions from m Q (H u ), m (Q 3 ) and m (U 3 ) 
approximately cancel out. Therefore, it is rather clear that splitting between these three scalar 
masses may substantially change the overall picture. If m (H u ) > ^/O.Qm (Q 3 ) 2 + 0Am (U 3 ) 2 
at the GUT scale, the scalars give negative contribution to fi 2 which makes the upper bound 
on mo(l,2) from the REWSB constraint more stringent unless gauginos are heavier. In the 
opposite case, the inverted hierarchy and the maximal stop mixing may be realized for lighter 
gauginos than in the case of degeneracy between mo(H u ) and universal third generation scalar 
masses. 



4 More phenomenology 
4.1 Dark matter 

In the IH scenario there are two distinctive ways to make the relic abundance of the LSF0 
compatible with the cosmological observations: 

• Stop-coannihilation for values of mo(l, 2) close to the boundary of the region where stops 
are tachyonic. The region of stop-coannihilation is generically present in the IH scenario 
for some intermediate range of mo(l,2) where the splitting between the stop NLSP and 
the neutralino LSP masses is small, ^dm^ 2 consistent with the WMAP bound is obtained 
for very small range of mo(l,2) and typically requires an adjustment of mo(l,2) with a 
precision of order 10~ 3 . This is because bino is the LSP which usually leads to too 
large values of ^dm/i 2 [35]. It is interesting to note that the stop-coannihilation region 
often coincides with the region where the stop-mixing correction to the Higgs mass is 
maximized. A benchmark point illustrating such a case is given as a Point A in Table [T] 

9 We use MicrOMEGAs [3J] to compute the relic abundance of the LSP, as well as BR(6 ->■ sj), BR(B S -> 
and SUSY contribution to muon anomalous magnetic moment, a^ USY . 
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Point A 


Point B 


Point C 


Point D 


Mi /2 


1000 


1500 


1500 


1000 


m (3) 


3700 


3400 


3800 


3300 


m (l,2) 


17690 


21070 


22500 


14500 




-2000 








-2000 


m (H d )/m (3) 


1 


1 


1 


1.6 


tan/3 


10 


10 


10 


50 


// 


888 


698 


452 


457 




125 


125 


125.1 


125.3 


m H 


3541 


3154 


3477 


3487 


m A 


3541 


3154 


3477 


3487 


m H ± 


3542 


3155 


3478 


3488 


A-1,2 
A3, 4 

m-± 


444, 813 


647, 707 


448, 461 


419, 467 


891, 940 


722, 1284 


677, 1286 


483, 869 


812, 940 


700, 1284 


455, 1286 


457, 869 


m 3 


2465 


3530 


3545 


2432 




17675, 17675 


21116, 21119 


22526, 22532 


14531, 14510 


™*1,2 


476, 1801 


699, 1581 


505, 1632 


979, 1274 


m d L ,R 

m b, 2 


17675, 17685 


21116, 21120 


22526, 22533 


14531, 14541 


1784, 2926 


1555, 2717 


1610, 2933 


1176, 1584 




17680 


21068 


22495 


14481 




3466 


3107 


3480 


2367 


m ~e L ,B 


17681, 17686 


21069, 21068 


22495, 22497 


14482, 14528 


m fl2 


3467, 3580 


3108, 3257 


3481, 3645 


1853, 2368 


&DMh 2 


0.111 


0.118 


0.021 


0.116 


BR(b ->■ s 7 ) 


2.92 x 10~ 4 


2.89 x 10~ 4 


2.66 x lO" 4 


1 x lO" 4 


BR(£ S ->■ //+//") 


3.07 x 10~ 9 


3.07 x lO- 9 


3.07 x 10~ 9 


3.61 x 10~ 9 


„SUSY 


1 x 10~ 12 


7 x lO" 13 


5 x 10~ 13 


1 x lO" 11 



Table 1: Several benchmark points with the inverted scalar mass hierarchy characterized by 
large stop-mixing contribution to the Higgs mass. Point A is an example of small mass splitting 
between the bino LSP and the lighter stop. Points B and D have a mixed higgsino-bino LSP 
Point C is characterized by mainly Higgsino LSP. 

• Higgsino LSP or mixed Higgsino-bino LSP for values of m (l, 2) close to the boundary of 
the region where /z 2 < 0. This resembles the focus point scenario in the CMSSM where 
for appropriately large values of m , |//| can become smaller than M 1 . In the present case 
small values of are obtained not only due to large values of mo (3) but also because of 
even larger values of mo(l, 2). If the Higgsino is the LSP its relic abundance usually turns 
out to be too small but if the LSP is a mixed state of bino and Higgsino the 2a WMAP 
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bound can be accommodated [35]. The region with a significant component of Higgsino 
in the LSP is much larger than the stop-coannihilation region and does not require very 
precise choice of parameters. However, the large stop mixing correction to the Higgs mass 
in the region with the Higgsino LSP is present only in some part of it where the Higgsino 
and the lighter stop are light simultaneously. In our plots the Higgsino and the lighter 
stop are both light in the part of the dark green region in the vicinity of the border 
between the grey and the yellow regions (corresponding to ji 2 < and tachyonic stop, 
respectively). The corresponding benchmark points are given as Points B and D in Table 
[I] for tan/3 = 10 and tan/3 = 50, respectively. 

4.2 b ->> S7 

The SM prediction [36] for BR(6 — > 57) is about la below the experimental central value |37j : 

BR SM (6 -»■ si) = (3.15 ± 0.23) x 10" 4 , 

BR exp (6 -»■ S7) = (3.55 ± 0.24 ± 0.09) x 10" 4 . (16) 

In the MSSM there are two important contributions to BR(5 — > 57) in addition to the SM ones: 
the charged Higgs contribution and the chargino-squark contribution [38J. The former one 
always increases the SM model prediction. In the IH scenario for tan (3 = 10 the charged Higgs 
is far above TeV so this contribution is negligible. On the other hand, the chargino contribution 
is negative (with respect to the SM) for > because the sign of this contribution is the 
same as sgn (fiA t ) [j^jand maximal stop mixing in IH scenario can be obtained only for A t < 0. 
It is also important to add that the main chargino contribution is maximized for the maximal 
stop mixing and grows linearly with tan/3. In the region of parameter space with the maximal 
correction to the Higgs mass from stop mixing BR(6 — > 57) is typically between la and 3a below 
the experimental central value for tan /3 = 10. The largest deviation from the experimental 
central value occurs for the points with the lightest Higgsino. 

Since the chargino contribution is proportional to tan /3, BR(6 — > 57) is even more restrictive 
for tan /3 = 50. In this case BR(6 — > 57) in the region of parameter space with the "maximal 
mixing" is far below the experimental value. For the benchmark point D in Table [T] with 
the Higgsino-bino LSP the discrepancy between the theory and experiment is about 8<r. This 
discrepancy can be relaxed to some extent if the Higgsino is heavier but even for ^ 1 TeV, 
BR(6 — > 57) is about 5a below the experimental value. 

4.3 B s ->> 

The LHCb upper limit for BR(B S ->• is [30]: 

BR(B S -> < 4.5 x 10~ 9 at 95%C.L. , (17) 

10 There is also part of the chargino contribution which has the same sign as sgn(— 11M2) so strictly speaking 
the chargino contribution may be negative also in some part of the parameter space where [iA t > 0. A detailed 
discussion on the sign of the chargino contribution can be found e.g. in [39] . 
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which is now very close to the SM prediction [41 J: 



BR SM (E S = (3.2 ± 0.2) x 10~ 9 . (18) 

This leaves very little room for contributions from the new physics. In the MSSM, BR(i? s — > 
A t+ A i_ ) probes the region of large tan/3 since the dominant MSSM contribution is proportional 
to A t tan 6 j3 /m\ [42]. At large tan/3, the CP-odd Higgs mass is given by: 

m 2 A ^m 2 Hd -m 2 Hu -M 2 z , (19) 

where m Hd and m Hu should be understood as the soft masses at the low scale. The RG running 
of m 2 Hd — rn 2 Hu at one-loop level is given by: 

8vr 2 ^(m| d -m 2 H J = 3h 2 X t - 3h\X h - h 2 T X T + p 2 S , (20) 

where 

X t = m 2 Q3 +m 2 U3 +m 2 Hu +A 2 , (21) 

X b = m 2 Q3 + m 2 D3 + m 2 Hd + A 2 b , (22) 

X T = m\ 3 + m| 3 + m 2 Hd + A 2 , (23) 



S is given by Eq. (12) and t = ln(M in /Q). We omitted terms proportional to the first and 
second generation Yukawa couplings which are negligible. Since the bottom and tau Yukawa 
couplings give the negative contribution to the RG running, the pseudoscalar Higgs becomes 
very light or even tachyonic when the top, bottom and tau Yukawa couplings are of the same 
order. In particular for tan/3 = 50 (and the other input parameters as specified below Eq. (|2])) 
m H d ~ 771 h u a t t ne scale Q — 1.5 TeV is given by: 

m 2 Hd - m 2 Hu w 0.6M 2 /2 - 0.02A M 1/2 + 0.03A 2 + 0.2m (3) 2 - 0.005m (l, 2) 2 . (24) 

The corresponding Yukawa couplings at the GUT scale equal: h t = 0.56, h b = 0.35, h T = 0.51. 
From the above formula it is clear that the CP-odd Higgs mass is driven to smaller values when 
mo(l, 2) increases. 

BK(B S — > /U + /i~) can easily be brought to phenomenologically acceptable values by splitting 
the soft Higgs masses at the GUT scale in such a way that m (H d ) > m (H u ) since such splitting 
makes the pseudoscalar Higgs heavier. 



4.4 MSSM spectrum and LHC phenomenology 

In the IH scenario with "the maximal mixing" the lighter stop is expected to be the lightest 
colored sparticle, in the range from 500 GeV to 1 TeV. The magnitude of the two loop effect of 
the heavy first two generations on the stau masses is controlled by the weak gauge coupling so 
they are typically heavier than the stops. Masses of the sfermions of the first and second gen- 
erations are set in the first approximation by mo(l, 2) so they are much heavier than sfermions 
of the third generation. 
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Gluino masses are in the range 2 — 3 TeV and other gauginos are correspondingly lighter, 
for universal M\/2- As explained earlier, higgsino mass may be in the range 200-400 GeV. 

If gluino is much heavier than the lighter stop, the stop pair production dominates SUSY 
production cross-section at the LHC. In such a case the limits on the stop mass are much weaker 
because the stop pair production cross-section is several orders of magnitude smaller than that 
of the gluinos with the same mass. The LHC limits on the direct stop production have not 
been presented so far. It is argued in jlEj that the existing LHC searches with jets and missing 
energy may have already excluded mi < 300 GeV in some particular cases when the Higgsino is 
the LSP and a main decay mode is t — > b^. However, this limit is highly model dependent and 
e.g. no lower mass limit for right-handed stop was found in [18] if bino is the LSP. Moreover, 
as stated before no official limits on direct stop pair production from the LHC experiments are 
available yetp] 

For large tan /3 the spectrum of the third generation is more compressed. In such a case the 
stop mass splitting is smaller so the lighter stop mass is expected to be a bit larger than for 
moderate tan (3 - not much below 1 TeV (in order to get Msusy ~ 1 TeV as preferred by the 125 
GeV Higgs mass), see benchmark point D in Table [l] In addition, sbottoms are expected to be 
in the TeV range. Constraining direct production of sbottoms is experimentally less challenging 
than that of stops so in IH scenario at large tan (3 sbottoms may be discovered (or ruled out) 
before the stops. Nevertheless, for the time being sbottom masses above about 400 GeV are 
consistent with the experiment [44] . 

From the above discussion it should be clear that the IH scenario is very weakly constrained 
at the moment. 

5 Conclusions 

The idea that the first two generations of sfermions are much heavier than the third one has been 
promoted in the past as a way to ease the supersymmetric FCNC problem, without violating 
the naturalness principle. Also, the sfermion mass spectrum can then be linked to the fermion 
masses in models based on horizontal symmetries. The early LHC results, putting stronger 
lower limits on the masses of the first-two-generation squarks than on stops and sbottoms, add 
some attractiveness to this idea. In this paper we have shown that this scenario predicts large 
stop mixing as a consequence of the RG evolution, with vanishing or small A-terms at the high 
scale. For the lightest stop mass to be at least C(0.5) TeV and assuming proper REWSB the 
Higgs boson is necessarily heavy, easily in the range 122-127 GeV. In particular, if the masses of 
the first-two-generation sfermions are above 15 TeV, the above conditions place a lower bound 
on the Higgs boson mass of about 122 GeV and the bound becomes more stringent as the masses 
of the first-two-generation sfermions increase. The Higgs boson mass of 125 GeV requires the 
universal gaugino mass to be about 1.5 TeV for vanishing A-terms at the GUT scale. It can 
be substantially smaller if negative A is assumed, e.g. for A = —2 TeV it is enough to have 
Mi/2 ~ 1 TeV. This scenario is only moderately fine-tuned. The LSP remains an interesting 

11 In |43j the ATLAS constraints on the direct stop pair production have been presented in a particular 
simplified model (motivated by gauge mediated SUSY breaking) with a gravitino LSP but the stop decay chain 
in such a model differs very much from those typical for the IH scenario. 
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dark matter candidate, particularly when it has a strong higgsino component or the stop NLSP 
is degenerate with the LSP to the extent which allows for efficient stop-coannihilations. The 
parameter range considered in this paper looks like a good bet for the MSSM. 



Acknowledgments 

This work has been partially supported by STFC. We would like to thank G. Gersdorff and 
R. Ziegler for useful discussions on inverted hierarchy models. MB would like to thank B. C. 
Allanach for helpful discussions, as well as T. Hahn, S. Heinemeyer and P. Slavich for useful 
correspondence. The present research was supported in part by the ERC Advanced Investigator 
Grant no. 226371 "Mass Hierarchy and Particle Physics at the TeV Scale" (MassTeV), by 
the contract PITN-GA-2009-237920 and by the French ANR TAPDMS ANR-09-JCJC-0146. 
S.P. acknowledges partial support by the contract PITN-GA-2009-237920 UNILHC and by the 
National Science Centre in Poland under research grant DEC-2011/01/M/ST2/02466. 



References 

[1] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 710 (2012) 49 [arXiv: 1202. 1408 [hep- 
ex]]; ATLAS-CONF-2012-019 

[2] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 710 (2012) 26 [arXiv:1202.1488 
[hep-ex]]; CMS-PAS-HIG-12-008 

[3] H. Baer, V. Barger, P. Huang and X. Tata, JHEP 1205 (2012) 109 [arXiv:1203.5539 
[hep-ph]]. 

[4] H. E. Haber, R. Hempfling and A. H. Hoang, Z. Phys. C 75 (1997) 539 [hep-ph/9609331]. 

[5] B. C. Allanach, A. Djouadi, J. L. Kneur, W. Porod and P. Slavich, JHEP 0409 (2004) 044 
[arXiv:hep-ph/0406166]. 

[6] P. Kant, R. V. Harlander, L. Mihaila and M. Steinhauser, JHEP 1008, 104 (2010) 
[arXiv: 1005.5709 [hep-ph]]. 

[7] P. Kant, arXiv: 11 11.7213 [hep-ph]. 

[8] B. C. Allanach, Comput. Phys. Commun. 143 (2002) 305-331. [hep-ph/ 01 04 145]. 

[9] G. Degrassi, P. Slavich and F. Zwirner, Nucl. Phys. B 611 (2001) 403 [hep-ph/0105096]; 
A. Brignole, G. Degrassi, P. Slavich and F. Zwirner, Nucl. Phys. B 631 (2002) 195 [hep- 
ph/0112177]; A. Dedes and P. Slavich, Nucl. Phys. B 657 (2003) 333 [hep-ph/0212132]; 
A. Dedes, G. Degrassi and P. Slavich, Nucl. Phys. B 672 (2003) 144 [hep-ph/0305127]. 

[10] M. S. Carena, H. E. Haber, S. Heinemeyer, W. Hollik, C. E. M. Wagner and G. Weiglein, 
Nucl. Phys. B 580 (2000) 29 [hep-ph/0001002]. 



21 



[11] M. Carena, S. Gori, N. R. Shah and C. E. M. Wagner, JHEP 1203 (2012) 014 
[arXiv: 11 12.3336 [hep-ph]]. 

[12] M. S. Carena, M. Olechowski, S. Pokorski and C. E. M. Wagner, Nucl. Phys. B 419 (1994) 
213 [hep-ph/9311222]. 

[13] M. S. Carena, M. Olechowski, S. Pokorski, C. E. M. Wagner, Nucl. Phys. B426 (1994) 
269-300. [hep-ph/9402253]. 

[14] N. Arkani-Hamed and H. Murayama, Phys. Rev. D 56 (1997) 6733 [hep-ph/9703259]. 

[15] F. Brummer, S. Kraml and S. Kulkarni, arXiv: 1204.5977 [hep-ph]. 

[16] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 192 (2011) 18 
[arXiv: 1001.4538 [astro-ph.CO]]. 

[17] Z. Berezhiani, P. H. Chankowski, A. Falkowski and S. Pokorski, Phys. Rev. Lett. 96 (2006) 
031801 [hep-ph/0509311]. 

[18] M. Papucci, J. T. Ruderman and A. Weiler, arXiv:1110.6926 [hep-ph]. 

[19] J. L. Feng, K. T. Matchev and T. Moroi, Phys. Rev. D 61 (2000) 075005 [hep-ph/9909334]. 

[20] K. L. Chan, U. Chattopadhyay and P. Nath, Phys. Rev. D 58 (1998) 096004 [hep- 
ph/9710473]. 

[21] S. Schael et al. [ALEPH Collaboration and DELPHI Collaboration and L3 Collaboration 
and ], Eur. Phys. J. C 47 (2006) 547 [arXiv:hep-ex/0602042]. 

[22] J. L. Feng, K. T. Matchev and D. Sanford, Phys. Rev. D 85 (2012) 075007 [arXiv: 11 12.3021 
[hep-ph]]. 

[23] S. Akula, B. Altunkaynak, D. Feldman, P. Nath and G. Peim, Phys. Rev. D 85 (2012) 
075001 [arXiv:1112.3645 [hep-ph]]. 

[24] F. E. Paige, S. D. Protopopescu, H. Baer and X. Tata, hep-ph/0312045. 

[25] M. Dine, R. G. Leigh and A. Kagan, Phys. Rev. D 48 (1993) 4269 [hep-ph/9304299]; 
A. Pomarol and D. Tommasini, Nucl. Phys. B 466 (1996) 3 [hep-ph/9507462]; R. Barbieri, 
G. R. Dvali and L. J. Hall, Phys. Lett. B 377 (1996) 76 [hep-ph/9512388]. 

[26] A. G. Cohen, D. B. Kaplan and A. E. Nelson, Phys. Lett. B 388 (1996) 588 [hep- 
ph/9607394]. 

[27] E. Dudas, S. Pokorski and C. A. Savoy, Phys. Lett. B 369 (1996) 255 [hep-ph/9509410]. 
E. Dudas, C. Grojean, S. Pokorski and C. A. Savoy, Nucl. Phys. B 481 (1996) 85 [hep- 
ph/9606383]; 

[28] C. D. Froggatt and H. B. Nielsen, Nucl. Phys. B 147 (1979) 277. 

22 



[29] P. H. Chankowski, K. Kowalska, S. Lavignac and S. Pokorski, hep-ph/0507133. 

[30] P. Binetruy and E. Dudas, Phys. Lett. B 389 (1996) 503 [hep-th/9607172]; G. R. Dvali 
and A. Pomarol, Phys. Rev. Lett. 77 (1996) 3728 [hep-ph/9607383]. 

[31] M. Gabella, T. Gherghetta and J. Giedt, Phys. Rev. D 76 (2007) 055001 [arXiv:0704.3571 
[hep-ph]]; O. Aharony, L. Berdichevsky, M. Berkooz, Y. Hochberg and D. Robles-Llana, 
Phys. Rev. D 81 (2010) 085006 [arXiv:1001.0637 [hep-ph]]; C. Brust, A. Katz, S. Lawrence 
and R. Sundrum, JHEP 1203 (2012) 103 [arXiv: 11 10.6670 [hep-ph]]; N. Craig, M. McCul- 
lough and J. Thaler, arXiv: 1203. 1622 [hep-ph], N. Craig, S. Dimopoulos and T. Gherghetta, 
JHEP 1204 (2012) 116 [arXiv:1203.0572 [hep-ph]]. S. Krippendorf, H. P. Nilles, M. Ratz 
and M. W. Winkler, Phys. Lett. B 712 (2012) 87 [arXiv: 1201.4857 [hep-ph]]. 

[32] S. P. Martin, In *Kane, G.L. (ed.): Perspectives on supersymmetry II* 1-153 [hep- 
ph/9709356]. 

[33] A. E. Nelson and D. Wright, Phys. Rev. D 56 (1997) 1598 [hcp-ph/9702359]. 
P. H. Chankowski, K. Kowalska, S. Lavignac and S. Pokorski, Phys. Rev. D 71 (2005) 
055004 [hep-ph/0501071]. 

[34] G. Belanger, F. Boudjema, A. Pukhov, A. Semenov, Comput. Phys. Commun. 176 (2007) 
367-382. [hep-ph/0607059]. 

[35] N. Arkani-Hamed, A. Delgado, G. F. Giudice, Nucl. Phys. B741 (2006) 108-130. [hep- 
ph/0601041]. 

[36] M. Misiak, H. M. Asatrian, K. Bieri, M. Czakon, A. Czarnecki, T. Ewerth, A. Ferroglia, 
P. Gambino et al., Phys. Rev. Lett. 98 (2007) 022002. [hep-ph/0609232]. 

[37] D. Asner et al. [Heavy Flavor Averaging Group], arXiv: 1010. 1589 [hep-ex]. 

[38] G. Degrassi, P. Gambino and G. F. Giudice, JHEP 0012 (2000) 009 [arXiv:hep- 
ph/0009337]. 

[39] M. Badziak, M. Olechowski and S. Pokorski, JHEP 1108 (2011) 147 [arXiv: 1107.2764 
[hep-ph]]. 

[40] R. Aaij et al. [LHCb Collaboration], arXiv:1203.4493 [hep-ex]. 

[41] A. J. Buras, Acta Phys. Polon. B 41 (2010) 2487 [arXiv:1012.1447 [hep-ph]]. 

[42] S. R. Choudhury and N. Gaur, Phys. Lett. B 451 (1999) 86 [hep-ph/98 10307]; K. S. Babu 
and C. F. Kolda, Phys. Rev. Lett. 84 (2000) 228 [hep-ph/9909476]. 

[43] ATLAS-CONF-2012-036 

[44] G. Aad et al. [ATLAS Collaboration], Phys. Rev. Lett. 108 (2012) 181802 [arXiv:1112.3832 
[hep-ex]]. 



23 



